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Abstract 

By a modular combination of mesoscopic detectors and gates, we present a thought- 
ful pump that transports gas particles against the difference of their density at the 
expense of the diffusion of another species of gas particles. We also discuss briefly 
the relevance of the model to the study of structure-function relationship of the 
biomolecular machines. 
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1 Introduction 



Biological molecular machines such as ion pumps, molecular motors (myosin, 
Fi motor, etc.) or signal transducers fG-p roteins) have a ttracted a lot of inter- 
est among physicists as well as biologists ({Alberts el al. . 1298). The molecular 



machin es works under th ermal fluctuations, where the latt er serve as an energy 
sourc e ( Sekimot'ol . 1997 ) for the free energy tra nsduction fe isen berg and Hilll . 



1981)2 as well as for the thermal activation (jWong et al.l . |1997l). Moreover 



the cycle of free-energy transduction in a single molecular machine requires no 
external control. As many of these molecular machines share highly conserved 
three-dimensional structural modules, it is widely believed that they have been 
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evolved from a common ancestors ( KuU et al. . 1998[ l. This fact suggests that 
such molecular machines may be composed of a few modular structures each 
of which bears some elementary functional role. With recent rapid develop- 
ment of nano-biology and the structural biology, the biologists have begun to 
look for such modular structure-function relationships within a single molec- 
ular motors ( Higuchi and Endow . 20021 : Nitta et al. . 2004 Tsiavaliaris et al. 
2004) . 



From the theoretical physicists side, there is no surprise about the mere 
fact that a molecular machine without any special forward-backward sym- 
metry can move in one direction under non-equilibrium con ditions (e.g. u n- 
der the ATP hydrolysis), as stated in the Curie's principle ( Curid . 1894 ) ^ . 
With this in mind, the question of our interest is how we can construct, 
in combining well-defined simple functional modules, an autonomously con- 
tro lled system that can wo rk un der thermal fluctuatio ns. Since the early works 
by Fevnman et al. ( 1966| ) and Biittiker et al. (1987), several models of au- 
tonomous free energy transducer have been proposed, and th eir improvements 
through the incorporat i on of a "gate" have been pr oposed (ISakaguchil . Il998l : 
Derenvi and Astumianl . Il999f ) (see also the review bv lReimannI (l2002[ )). These 
models were, however, rather ad hoc and the modular nature of the construc- 
tion was not clearly visible. In the present paper, we would construct a free 
energy transducer with being more conscious about the modules and taking 
the object oriented approach (|Sekimotol . l2000h . We start by introducing the 
concept of "semi-detectors", the detectors that can perceive a certain exter- 
nal state, but not all, with arbitrary sureness despite the thermal fluctuation 
(§ 2). We then combine these semi-detectors together with the gates, under 
the designing concept that we call the "bidirectional control", to realize an 
autonomous particle transporter (§ 3). Because of the limited space of this 
special issue, the extensive discussion on the possible relations to the real 
biological molecular machines will be given elsewhere. 



2 semi-detectors in the fluctuating world 



The detector is the mechanism that correlates the states of the outside to those 
of the inside of the system in question. We focus on the case where the de- 
tection site admits a single ligand particle coming randomly from the outside. 
We will limit ourselves to the energetically "neutral" detections, involving no 
change of the total energy. Let us define the mapping from a state of the 
outside, X, to a subset of the states of the inside, $(x). In the discrete repre- 
sentation, the states of outside consists of IN (i.e. a ligand particle is on the 

^ "Lorsque certaines causes produisent certains effets, les elements de symetrie des 
causes doivent se retrouver dans les effets produits. " 
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detection site) and OUT (otherwise), while those states inside consist of ON 
and OFF. The perfect detector would establish the mapping, $(IN) = {ON} 
and $(OUT) = {OFF}. But these correspondences are too stringent to realize 
under thermal fluctuations. We will show below that there can be the physical 
mechanisms which assure either one of the above correspondences. We will 
call the mappings corresponding to such mechanisms the semi-detectors. 

Semi- detector of absence: This module functions as the mapping $ab(2;) that 
prohibits only the output of OFF under the input of IN. That is, 

$ab(IN) = {ON}, $ab(OUT) = {ON, OFF}. (1) 



Then the OFF surely indicates the absence of the particle (OUT). Such aspect 
is useful for r epressive pro cesses like the suppression of DNA transcription by 
an repressor ( Lewin . 2003| ). In the continuum representation, one may define 
X as the position of a ligand particle in a half space <x< oo, where a; = 
corresponds to the detection site (see Fig. 1 Left). The state of the inside, a, 
is assumed to be bounded, say, in the region of — 1 < a < without loosing 
generality. The semi-detector of absence may be realized by the steric repulsion 
between the ligand particle and a movable object (thick bars in Fig. 1 Left). 
The total energy of the ligand-detector system, f/(x, a), is written as ?7(x, a) = 
for x> a + l and U{x,a) = oo otherwise. OFF and ON states correspond, 
respectively, to a>—l + 6 and — l<a<— 1 + 5 with a small 6>0. 

OUT 



ON 



ON 



OFF 
IN 








+1 



Fig. 1. Left: Schematic representation of the semi-detector of absence in a continuum 
representation. The movable object (thick bar) and the ligand particle (thick dot) 
are sterically excluding with each other, a represents the coordinate of the left 
extremity of the object, x = is the detection cite. Right: The shadowed region on 
the (x,a)-plane indicates the accessible phase space, where U{x,a) =0. 



Semi- detector of presence: This module functions complementarily to the for- 
mer semi-detector. The mapping $pr(a;) that it defines prohibits only the out- 
put of ON under the input of OUT. That is, 

$pr(IN) = {ON, OFF}, <l'pr(OUT) = {OFF}. (2) 



Then the ON surely indicates the presence of the particle (IN). Such aspect is 
useful for active processes like the uptake of ATP by a molecular motor. In 
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the continuum representation, the semi-detector of presence may be reahzed 
by the compensation of a strong restoring potential for the movable object, 
M0(a) (M^^/cbT), by the strong attractive interaction energy, — M0(a — x), 
between the ligand particle and the movable object (see Fig. 2 Left). Here (f){z) 
is defined on — oo <z<0 so that (j){z) = for 2;<— 1 and that (f){z) increases 
monotonically from </>(— 1) = to 0(0) = 1. Energetically, this movable object 
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Fig. 2. Left: Schematic representation of the semi-detector of presence, a and x 
represent, respectively, the position of the movable object (filled square) and that 
of the ligand particle (thick dot). Unless the particle is within the proximity of the 
detection cite {0 < x < 6) with small 6, the movable object is constrained at a ~ — 1 
even under thermal agitations. Right: The accessible phase space is indicated as a 
shadowed region, where U{x, a) ~ k^T. 

can take the value a > —1+6 only when the particle is near the binding site, 
that is < X < 5 with a smaU 5 > so that M(f){6)<kBT^ . Therefore ON 
and OFF correspond, respectively, to a > — 1 + 5 and — 1 < a < — 1 + 5. The 
compensation mechanisn a similar to the o ne discussed here has been discussed 
and called induced fit by Koshland ( 1973[ l in the context of the ligand binding 



by enzymes. 

The above two semi-detectors establishes the correlations between the outer 
world of the detectors and the detectors themselves in the way that the 
semi-detectors represent certain error-free information despite the thermal 
fluctuations. If the state variables of two semi-detectors of presence, say Oi 
and 02, respectively, are coupled energeticaUy with each oth er, the result- 
ing system might function t o realize the "cooperative bin ding" (iMonod et al. 



1965 ; Koshland et al. , 1996t ) or the "exchange of binding" ( Eisenberg and Hill 



19851), depending on the nature of the coupling. 



3 Free-energy transducer based on the semi-detectors of presence 
and the gates 



In order to focus on the subject of the autonomous control, we would like to 
avoide energetic aspects as far as possible. To this end, we will construct a 
pump of the ideal "load"(L) gas particles from a dilute reservoir, (L,l), to a 

^ One can verify this by drawing M{(j)(a) — (j){a — x)} vs a for various values of x. 
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dense reservoir, (L,h), at the expense of the diffusion of the ideal "fuel"(F) 
particles from a dense reservoir, (F,h), to a dilute reservoir, (F,l), see Fig 3. 
If we denote the chemical potentials of these reservoirs by /iL,h /^L,h, /^F,h and 




Fig. 3. Schematic diagram of how the free-energy transducer (the hexagon) pumps 
up the load particles at the expense of the diffusion of the fuel particles. See the 
text. 

yUp,!, respectively, then the decrease of the total Gibbs' free energy in the course 
of the transport of AiV^ of the load particles at the expense of ANp^^^i of 
the fuel particles is (/xp.h — /^F,i)A^F,h^i— (/^L,h — /^L,i)A^L,wh- Below we will 
present a completely symmetric model in which the roles of the fuel and the 
load are totally exchangeable depending on the relative magnitude of /xp.h— /^f,! 
and /XL,h - /iL,i ^ • 

The cyclic process in the above mentioned pump requires at least two internal 
degrees of freedom ^ . To construct such pump with using the semi-detectors 
of presence as constituent modules, the most natural way would be to put one 
semi-detector on the fuel side (with the variable ap) and the another one on 
the load side (with the variable cl), and to have them to con trol the access 



of th e particles on the opposite sides, via allosteric couplings (jMonod et al 
see Fig. 4. We would call this scheme "bidirectional control". We will 



describe in three steps the details of the model which implements this scheme: 
Step-1: Introduction of the "reaction" coordinates Xl andxp. Being inspired 
by Fig. 2, we introduce the convenient coordinate x (0 < a; < oo) which 



^ Under this symmetry, the currents of the active transport of the load particles, 
JL,i^h) and that of the passive diffusion of the fuel particles, Jp.h^h ai'e given in 



A/Usym- Here Lasy > and Lsym ^ are the kinetic 



the linear non-equilibrium thermodynamics as Jl,i 

and JF,h-+l = -^^asy^/^asy + ^ 

coefficients, and A^asy = AtF,h-/^F,i-A*L,h+/iL,i and A/isym = MF,h-MF,i+AiL,h-ML,i- 
This pump can work almost reversibly near the stalled condition, A/ig^y = 0, if Lgym 
is sufficiently small. 

^ This state ment does not contradicts with the f:act tha t the existing models of 
heat engine ( Fevnman et al. . 19661 : Biittiker et al. . Il98?f ) have assumed only one 
rotational degree of freedom to lift a load, since such degree of freedom corresponds 
to two bounded degrees of freedom. 
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Fig. 4. Left: Scheme of the bidirectional control. Right: Typical cycle realized by the 
two independent internal state variables. A-F correspond to those in Fig. 7. 

describes a ligand particle and a semi-detector of presence: x = max(5— 1, 0), 
a = — min(x, 1). Then we extend this definition of x to the case of two particle 
reservoirs: We assume that the high-density [low-density] reservoir occupies 
the half space x > [x < 0], respectively. We redefine the mapping x i— > {x, a) 
so that X can take the values on the entire axis, — oo < x < oo, with: 



X = {x/\x\) max(|a;| — 1, 0), a = — min(|a;|, 1). 



(3) 



A ligand particle is on the detection site if |x| < 1, in the high-density reser- 
voir if 5; > 1, and in the low-density reservoir if 5; < —1. We interpret x so 
that the region oi \x\ < 1 — 5 with small 5 > corresponds to the state ON 
of the semi-detector, while |x| > 1 — 5 corresponds to the state OFF. Finally 
wc apply this type of mapping for both the load side (xl) and the fuel side 
(xp). The coordinate plane {xi^,xp) can then represent the positioning of one 
representative load particle and the one fuel particle together with the states 
of the semi-detectors of presence, cl and ap ^ . 

Step-2: Definition of the gates' action. For the fuel particles, we construct 
the gate which allows the access of the particles exclusively from one of their 
reservoirs at a time. In the Top of Fig. 5, the potential barrier (the hight 
^ ksT) is established at xp = 1 + e with a small e > 0, so that the access 
of the fuel particles from the high-density reservoir is blocked at a distance e 
off the detection site {\xf\ < 1). The same architecture is defined for the load 
particles. Similarly the gate represented in the Bottom of Fig. 5 blocks the 
access of the fuel particles from the low-density reservoir at the distance e off 
the detection site. 

Step-3: Coupling of the semi-detectors to the gates. We define the bidirec- 
tional control by the following symmetric rules: 

[Control by a-p] If a fuel particle is detected, i.e. il\xp\ <1 ~ 5 {Top-Left of 
Fig. 6), only the load particles in the low-density reservoir (xl < ~1) can access 
to their detection site. If a fuel particle is not detected i.e., if |5;f| > 1 — 5 {Top- 

By the introduction of xl and xp^ we can avoid to be lost in the four-dimensional 
space. The trade-off that we pay is that any value of ap > — 1 is represented by the 
two equivalent points, xp = ztap. The same holds also on the load side. 
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Fig. 5. Top: Potential profile of the gate that opens to the fuel particles in their 
low-density reservoir. Bottom: Similar to the above but the accessibility is for the 
high-density side. The thick horizontal bars indicate the detection site for the fuel 
particles. 




Fig. 6. Top: The control of the load particles by the detection of the fuel particles. 
Bottom: The control of the fuel particles by the detection of the load particles. The 
J-shaped symbol in, for example, the Top-Left indicates the exclusive accessibility 
from the low-density reservoir of the load particles, as indicated in the Top of Fig. 5. 
The other cases would be understood similarly. 

Right), then only the load particles in the high-density reservoir (xl > — 1) 
can access to their detection site. 

[Control by ai]: If a load particle is detected, i.e. if < 1 — S {Bottom- 
Left of Fig. 6), only the fuel particles in the low-density reservoir {xp < — 1) 
can access to their detection site. If a load particle is not detected, i.e. if 
> 1 — S [B ottom- Right) , then only the fuel particles in the high-density 
reservoir {xp > — 1) can access to their detection site. 

The consequence of these simple and symmetric combination of the semi- 
detectors and the gates is immediately seen by the graphical representation 
in Fig. 7. There, the potential barriers of the gates are indicated by the thick 
horizontal or vertical bars. For example, the short horizontal bar at xp = 1 + e 
in the range of I^lI < 1 — S represents the blockade of the access of the 
fuel particles from their high-density reservoir. We see how the combination 
of those bars organizes a broad passageway joining the second and the forth 
quadrants on the (xl, xp)-p\a.ne. This passageway corresponds to the diffusion 
of a single fuel particle from its high-density reservoir to the low-density one 
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Fig. 7. In the (xl, XF)-plane, the thick horizontal and vertical bars represent the 
presence of the potential barriers either for the fuel particles (the horizontal bars) or 
for the load particles (the vertical bars). The shaded regions arc where the bidirec- 
tional control prohibits to access except for a small "leak" due to the finite values 
of e and 6. See the text for the details. The dashed curve shows one representative 
process (cf. the footnote on the multiplication of the state points). B-E on the curve 
correspond roughly to those in Fig 4. 



accompanying the active transport of a single load particle from its low-density 
reservoir to the high-density one. Therefore, if we could neglect the small gaps 
between the bars due to the finite values of e and S, this pump would work 
tightly with AA^^h-»i = ^^L,i-+h) whose mean flow direction is determined by 
the sign of (/xp.h — A*f,i)— (A*L,h — A'-l,!)- Actual autonomous pump seems not 
be able to avoid the "leak" due to the above mentioned gaps, but this leak 
would be of importance only near the stalled state (//F,h — A*f,i = A*L,h — A*l,i) 
as far as e and 5 are small ^ . 

How would the pump thus constructed looks like for an observer who can 
survey only the load particles? Let us suppose that (/iF,h — A*f,i)— (A*L,h — A*l,i) 
is positive and enough greater than UbT. As long as the detection site of 
the load particle is empty, this site is almost always accessible from the low- 
density reservoir of the load particle (B along the dashed curve in Fig. 7). 
When a load particle arrives at the detection site and is detected (C,D), the 
gate is very likely to reverse its accessibility so that the load particle on the 
detection site can leave now for the high-density reservoir (E). Then as soon 
as the load particle quits the detection site (F), the gate comes back to the 
initial conformation (B). Thus the pump behaves as if it responded by itself 
to the arrival and the leave of the load particle. 



° In the context of the the linear non-equilibrium thermodynamics, we expect that 
-^syiii/^asy IS at most of the order of e and S. 
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4 Discussion 



We have constructed theoretically an autonomous system that works under 
thermal fluctuations. The errors of the detection was avoided by introducing 
what we call the semi-detectors. The concept of the semi-detectors on meso- 
scopic scale might also be of interest in the context of the mesoscopic devices 
(jlmxiJ, |l997|) . The cyclic process was enabled by what we call the bidirectional 
control. The latter idea is applicable also to the macroscopic autonomous 
processes. For example, in the operation protocols of vending machines or 
of pay-phones, the gate on one side (taking up the money) and that on the 
other side (rendering goods or services in exchange) are controlled by their 
respective detectors on the opposite sides. 

The biological molecular motors are not as simple as we have discussed. Still, 
the present analysis might help as a reference frame when we look for the 
structure-function relationship in those systems. For example, we might ask 
if the hydrolysis reaction of ATP (fuel) in a molecular motor corresponds to 
changing the accessibility from the reservoir of the ATP to that of the ADP and 
the inorganic phos phate [Pi]. Also we might wonder if the putative Pi-sensors 
( "switch loop" etc.. Valj ( 1996t ): Brendza et al. (2Qq3)) can be compared to the 
semi-detector of the ligand in our model. As for the kinetics of the interaction 
between the molecular motor and its counterpart filament (load), we might 
ask if the motors have a degree of freedom to detect its own mechanical strain 
f)Tsiavaliaris et al.L 12004 : ISekimotd . liiooh in addition to the degree of freedom 
to control the binding to the filament (jNitta et al.l . 12004^ . As a qualitative 




"Wilde-type" 



"Mutant" 



Fig. 8. The pump of the "wild-type" (left, see Fig. 7) and a "mutant" with blocked 
fuel reactions (right). See the text. 

prediction, we might mention the possibility of the "mutants" in which the 
degree of freedom related to the detection of, say, the load particles, ol, is 
immobilized, as shown in Fig. 8. In this case, the rate of passive diffusion of 
the load particles should depend on the density of the (non-consumed) fuel 
particles in their high-density reservoirs. 
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